Complete Adiabatic Quantum Search in Unsorted Databases 
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We propose a new adiabatic algorithm for the unsorted database search problem. This algorithm 
saves two thirds of qubits than Grover's algorithm in realizations. Meanwhile, we analyze the time 
complexity of the algorithm by both perturbative method and numerical simulation. The results 
show it provides a better speedup than the previous adiabatic search algorithm. 
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Quantum computation is a promising way to solve clas- 
sical hard problems. Several quantum algorithms have 
been designed to perform classical algorithms with re- 
markable speedups. The most useful one among these 
is Grover's algorithm[l] which concerns the problem of 
searching for a required item in a unsorted database. 
One common example for this unsorted database search 
is to find a person's name in a phone book (the items are 
sorted by names) with only knowing his phone number. 
Classically, the only way to achieve this is brute-force 
search[l, 2] which requires an average of ^ quires for N 
entries in the phone book. However, if the information 
is stored in a quantum database, to find the right name 
with Grover's algorithm costs only a time of order \/iV, 
providing a quadratic speedup. 

The main process of Grover's algorithm is, swing- 
ing the index qubits from an initial uniform state to 
approach the solution state. The information of the 
database is not explicitly accessed in the processing of 
search. Instead, an Oracle is supposed to know all the in- 
formation in the database and act properly towards a in- 
put state depending on whether it denotes the solution [3]. 
Early experiments[4, 5, 6, 7, 8] of Grover's algorithm 
constructed this Oracle from a marked state as a func- 
tional analog instead of querying the database. For a 
complete solution of search problem, Kim and cowork- 
ers proposed a new approach to realize the Oracle on a 
quantum database and implemented it in experiment [9]. 
In this complete approach of Grover's algorithm, extra 
qubits are used to store the database. Similar method 
for construction of Grover's Oracle was also theoretically 
discussed later [2]. 

While Grover's algorithm is presented in the standard 
circuit model(i.e., using a sequence of discrete quantum 
gates), a new model of quantum computation shows up 
where the state of quantum computer evolves continu- 
ously and adiabatically under a certain time-dependent 
Hamiltonian. This new adiabatic model was soon ap- 
plied to the database search problem[l()] and the original 
adiabatic search algorithm was proved to have a time 
complexity of order N, which is the same performance as 
classical algorithms. More recently, Roland and Cerf[ll] 



recovered the advantage of adiabatic search to order v N 
(the same with Grover's algorithm), by performing the 
adiabatic evolution locally. However this adiabatic search 
algorithm construct the Hamiltonian from a marked state 
instead of referring to the database, thus it is not a com- 
plete search algorithm. And in order to discriminate it 
from our algorithm, we call it the marked-state adiabatic 
sefflrc/i(MSAS) algorithm in the following paragraphs. 

In this letter, we apply the quantum adiabatic com- 
putation to unsorted database search problem again and 
present a new quantum search algorithm. We will put 
forward a new method to represent the database. By 
this method the algorithm contains no Oracles and saves 
I of qubits than the complete approach of Grover's 
algorithm[' ]. We also analyze the time complexity by 
both perturbation method and numerical simulation. 
The results show it provide a higher speedup than the 
MSAS algorithm. 

As a new quantum computation model, adiabatic al- 
gorithm was brought out by Farhi et al.[12\ and soon 
became a rapidly expanding field. The idea of this new 
computation model is to prepare a system in the ground 
state of a simple initial Hamiltonian, then slowly switch 
the simple Hamiltonian to a complex Hamiltonian whose 
ground state encodes the solution to the problem of in- 
terest. According to Adiabatic Theorem, the the system 
stays in the ground state of the instantaneous Hamil- 
tonian if we perform the evolution slowly enough. So 
finally the state describes the solution to the problem. 
The time-dependent system Hamiltonian is 



H{t) = [1 - s{t)]H, + s{t)Hp, 



(1) 



where Hi is the initial Hamiltonian and Hp is the problem 
Hamiltonian which encodes the solution, and the mono- 
tonic function s{t) fulfills s(0) = and s(T) = 1. 

Here let's focus on the unsorted database search prob- 
lem. To be simplified, the database is a list of (i, Vi) pairs 
and sorted by i where i denotes index and Vi is value. 
Both i and Vi are rt-bit binary codes thus the database 
contained TV = 2" items. The "unsorted" property of 
the database refers to the field value not index. The 
unsorted database search problem here is looking for the 
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corresponding index i for a given target value t. And we 
assume that there's only one solution in the database for 
each search. Next we will describe the process to find the 
right i which connects to the target t. 

The essential part of an adiabatic search algorithm is 
how to encode the solution in the ground state of problem 
Hamiltonian. For example, the MSAS algorithm con- 
structs the problem Hamiltonian as Hp = 1 — |m)(m| 
where |m) is exactly the solution state. Thus it is 
not a complete database search. Obviously for a com- 
plete search, the information in database should be rep- 
resented in quantum forms. Taking the complete ap- 
proach of Grover's algorithm as an example['2, 9], the 
database is represented in an operator which satisfies 
= Uf generates the entanglement of 

qubits to denote the relation between i and Vi, thus both 
the fields are represented by qubits. 

In the present algorithm, however, not both the fields 
are represented by qubits. We define a database operator 
as 

JV-l 

V=J2vS){^\■ (2) 

i=0 

Clearly in this approach, the index is represented by 
qubits while the value is store in the strength of interac- 
tions. So no extra qubits are needed for the database. 
The operator 2? contains all the information in the 
database. Thus, we can construct the problem Hamil- 
tonian from V as 

Hp = {V~t)\ (3) 

where t is the target value which we are looking for. 

To test the validity of Hp, we will examine its ground 
state. To this end, we can write iJp as^^Q^(wi— 
From this form, each diagonal element of Hp is the square 
of difference between Vi and t. Thus the ground state will 
be the solution state \i) where Vi equals to t. 

Of course this construction provides a valid prob- 
lem Hamiltonian for the search problem. However, the 
Hamiltonian in Eq.(3) has a spectral width exponentially 
growing with the number of qubits, which is hard to re- 
alize when the database is large. Thus it is only useful 
in small-size databases. 

To solve this problem, we divide the comparison be- 
tween Vi and t into n sub-comparisons, each of which 
is performed for a single bit between them. Thus the 
database operator should be formed separately for each 
bit. For the jth bit of value we define the bit database 
operator "Dj as 

where vij is the j th bit of Vi . Similarly with the operation 
in Eq.(3), the problem Hamiltonian for each bit is 

Wp = (V,-tjf, (5) 



where tj is as well the jth bit of t. Consequently, the 
overall problem Hamiltonian is the summarization of all 
bit problem Hamiltonians 

n— 1 n—1 

Hp = J2H^ = J2[V,{l~t,) + t,{I-V,)] 

3=0 j=0 
n-1 

where tj is the complementation of binary bit tj and 
Also for a test of the validity, we can simplify Hp as 

Hp = (7) 

i=0 

where the function hivi^t) is the Hamming distance be- 
tween Vi and t. Thus the state \i) where h{vi,t) — 
is the ground state of Hp and is also the solution state. 
Moreover, the spectrum was successfully bounded in a 
range from to n. 

After the preparation of problem Hamiltonian, we will 
choose an initial Hamiltonian Hi . Hi should be chosen to 
be noncommutative with Hp to avoid crossing of energy 
levels[l()]. Normally, Hi is 

H, = gK + Tl + '-' + T^-'), (8) 

which means the qubits coupling with a magnetic field 
at the x-direction and the coupling strength is g. The 
ground state of Hi is 

where b{j) is the Hamming distance between j and 0. 

In the adiabatic evolution, the system Hamiltonian in- 
terpolates from Hi to Hp (i.e., see Eq 1) and the state 
of the system evolves according to the Schrodinger equa- 
tion. If this evolution acts adiabatically, the system will 
always stay on the instantaneous ground state of H{t) 
and in the end the solution of our problem will show up. 

An explicit application is necessary for a clear under- 
standing. Here we perform a 3-bit unsortcd database 
search for example. We randomly generate a database in 
a list as {6, 3, 5, 0, 4, 1, 7, 2}. The position of each value 
in the list refers to the index which ranges from to 7. 
For convenience we rewrite the database as binary codes 
which is {110,011,101,000,100,001,111,010}. Because the 
database operators and problem Hamiltonian are diago- 
nal, they are expressed by only the diagonal elements. 

X>o ^ dmg{0, 1,1, 0,0, 1,1,0} 
Vi ^ dm.g{l, 1,0, 0,0, 0,1,1} 
V2 = (img{l,0, 1,0, 1,0, 1,0}. (10) 
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FIG. 1; Process of the adiabatic evolution to search for 
value 5 in the mentioned database. (a)The instantaneous 
eigenvalues of the system Hamiltonian as a function of s. 
The solid hue represented the energy level of ground state, 
(b) Occupation probabilities of the system for the computa- 
tional basis during the adiabatic evolution in the numeri- 
cal simulation. The system starts from a uniform state and 
evolves to the solution state |010) which shows value 5 is on 
index 2. The parameter g of Hi in this example is 0.5. 



After building the quantum database, the problem 
Hamiltonian can be constructed for each search task. For 
example, if we want to find the position of the value 5 
which is 101 in binary, the problem Hamiltonian is 



Hp = V0 + V1+V2 

= rfm.g{2,2,0,2,l,l,l,3}. 



(11) 



To perform the adiabatic evolution, we initially pre- 
pare the system on the state of Eq.(9), and adiabatically 
switch the system Hamiltonian from a initial Hamilto- 
nian in Eq.(8) to the problem Hamiltonian in Eq.(ll). 
Finally the system will on the ground state of Hp which is 
the state |2). After measurement, we can get the knowl- 
edge that value 5 is on position 2. Fig.l shows the process 
of the adiabatic evolution for this search. 

For the practical usefulness of an algorithm, the occu- 
pied amount of resources is an important aspect. With- 
out doubts, the number of qubits needed for our algo- 
rithm equals to the bit width of the index because only 
the index field is represented by qubits. Thus the spa- 
tial complexity is n. As a comparison, since both the 
fields value and target are represented by qubits, the 
spatial complexity of the complete approach of Grover's 
algorithm [9] is 3n. Although the MSAS algorithm also 
have a spatial complexity of rt, it is not a complete 
database search. Thus our algorithm has the best spa- 
tial complexity in the quantum algorithms for complete 
database searches. 

To evaluate the time complexity of our algorithm, a 
decisive mathematical analysis is not possible. Therefore 



in this letter, we use both the perturbative method[13] 
and numerical simulation[12] to examine the situation of 
the time complexity. 

In the perturbative approach[13], the time cost of a 
adiabatic algorithm by either global or local evolution 
can be written as 



Tlocal « v/l^(X y/\S~\/\S+\, (12) 

log 1/(5 



S w {z : h{z, /) < mc\, nric oc 
S- w {z:E,<E,],E,^ 



logC^ 
logl/C" 



(13) 



where S'^ is a set containing the eigenstates of the prob- 
lem Hamiltonian which have a small Hamming distance 
towards the solution state |/), while S" contains the ones 
which have low energy levels. \S\ is the cardinality of set 
S . arc dimensionless parameters which are defined as 
C,^ = C(s* ±eo)- Here C(t) — i^lj^tj and s* is the position 
of the minimum gap between the ground and first exited 
state, eg and S are small numbers. 

To apply this result to our algorithm, we assume that 
the minimum gap is on the central position of s, thus we 
can get = 1/C~ > 1- Then we define a small number 
n = ■ Because the degeneracy of energy levels in 

the problem Hamiltonian in Eq.(7) is C*j where i is the 
ith energy level, Eq.(13) goes as 



i=0 

\S-\ ^ J2 Q,i?cCXl/0. 



(14) 



i=0 



Here, f2 is a small number and is not supposed to increase 
with 71, so only some low energy levels will be in 5"" 
and 15^1 is a comparatively small positive integer. Since 
J27=o ~ ^1 loi' 11^*3 worst case, we can take Eq.(12) 
as Tgiobai oc \S~\ oc N" where a < 1 is a constant. 

To derive a more accurate range for a, we per- 
formed a numerical simulation [12] for randomly gener- 
ated databases with the bit width of index sized from 
5 to 16. For each bit, we randomly generated 50 in- 
stances of database search. Then we performed a nu- 
merical global evolution using four-order self-adapted 
Runge-Kutta method to get a success probability of range 
[0.12, 0.13] for each instance. The mean time for each bit 
is shown in Fig. 2. By fitting the mean time, we obtain 
a = 0.81. For a comparison, we simulated the time com- 
plexity of MSAS algorithm using the same environment. 
The value a of MSAS algorithm obtained from the fitting 
is 1.02. The result of simulation fit well with theoretical 
expectation where a is 1[1()]. 

In Fig. 2, the running time of our algorithm grows 
much more slowly than the MSAS algorithm. This result 
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FIG. 2: Comparison of running time between our algorithm 
and MSAS algorithm to achieve a success probability of 1/8 
as the function of bit width of index. The solid line is the fit 
of the triangles each of which represent the mean simulation 
result of 50 instances (32 instances for n = 5) of our algorithm, 
while the dash line is the fit of the circles each of which is 
that of the MSAS algorithm. The error range of fitting the 
triangles is ±2% and that for the circles is ±0.3% 



matches well with the expectation from the perturbative 
analysis. Both the results show that our algorithm has 
a better performance in time complexity than the MSAS 
algorithm. And because local evolution can provide a 
quadratic speedup over global evolution, theoretically the 
time complexity of our algorithm by local evolution can 
be reduced to less than order \/N, even lower than than 
the complexity of Grover's algorithm. 

To be concluded, we introduce a new algorithm for 
quantum search problem by adiabatic evolution. Wc use 
another method to represent the quantum database in 
this algorithm and it saves | of qubits than the com- 
plete approach of Grover's algorithm [:)]. Wc use both 
the emerging perturbative method of adiabatic algorithm 
and numerical simulation to analyze the time complex- 



ity in this algorithm. The results show that it provides 
a higher speedup than the MSAS algorithm and poten- 
tially has a better performance than Grover's algorithm. 
This algorithm can be experimentally verified in NMR or 
ion-trap systems [14, 1-5]. 
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